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The hybridization of the single- excitation branch with the two- excitation con- 
tinuum is reconsidered from the theoretical point of view by including the 
effect of the interference term between one and two excitations. The phe- 
nomenological theory presented is used to reproduce the experimental data 
over a wide region of momentum (2.3-3.6 A' 1 ) and energy (0-12 meV). It 
is thus possible to extract the final part of 4 He spectrum with higher accuracy. 
It is found that data agree with a dispersion relation always below twice the 
roton energy. This is consistent with the negative value of the roton-roton 
interaction found. 

PACS numbers: 67.40. Db, 61.12.-q 

INTRODUCTION 

Excitations in superfluid 4 He have been widely studied in the last decadesJ 
However, the nature of the single particle spectrum termination is still un- 
clear. Pitaevskii a long time agcll predicted a termination of the spectrum 
due to a decay into pairs of rotons. Qualitatively the theory predicts that the 
low energy pole is repelled by the continuum, so that the spectrum flattens 
out for large Q towards 2 A losing spectral weight (A is the roton energy). 
At the same time, a damped excitation for u > 2A appears and shifts to 
higher energies. More precisely, the end of the spectrum is determined by 
the sign of the roton-roton (V4) interaction. If this is positive, the spec- 
trum terminates non analytically at 2 A and at a definite momentum Q c . 
On the contrary, if V4 < the spectrum should go smoothly to the two- 
roton bound state energy slightly below 2A. In this second case a sharp 
peak should survive up to large values of momentum, even if its spectral 
weight should decrease very rapidly. (This second scenario was proposed by 
Zawadowski-Ruvalds-SolanatH and Pitaevskii.a) In both cases no damping of 
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the excitations is possible below 2A, nor the sharp excitation can exceed 2A. 

Neutron scattering experiments suggest that the decay of excitations 
into pairs of rotons does actually take place for momentum Q > 2.6 A _1 .H 
Despite the good qualitative agreement between theory and experiment the 
sign of V4 is still undetermined and the possibility that the excitation energy 
exceeds 2A seems suggested by the data. Direct fitting of Pitaevskii-ZRS 
theory to data! indicates that V4 > 0, but experimental finding of a quasipar- 
ticle peak above 2A is not accounted for by the theory with reasonable values 
of the parameters. As a matter of fact, the position of the quasiparticle en- 
ergy was extracted by fitting a Gaussian peak on a background of constant 
slope.! Furthermore, the experimental finding of a (large) positive V4, from 
the theoretical point of view, is completely inconsistent with the observation 
of a spectrum that does not disappeaxat any finite Q. Recent experimental 
investigations by Fak and coworkers^ on the temperature dependence of 
the dynamical structure factor S(Q,tu) shows clearly that there is a strong 
correlation between the low-energy peak and the high energy continuum as 
Q increases from 2.3 to 3.6 A _1 .0 Thus supporting the hybridization picture. 

In this paper we address these inconsistencies that we believe are mainly 
due to the difficulty of taking properly into account the contribution of the 
continuum of excitations starting at 2A, when, for large Q, it becomes more 
important than the discrete contribution of the sharp state slightly below 
2A. 



1. PITAEVSKII-ZRS THEORIES AND THEIR EXTENSION 

We recall that the validity of Pitaevskii and ZRS theories is restricted 
to a small region around 2A. Indeed Pitaevskii in his original paper@ ex- 
ploited the logarithmic divergence appearing in the two-roton response func- 
tion [F (Q,u) = i /ff /|^S(„0£Cp'), where S^CP) =u+iO+-u(Q), 
w(Q) is the measured spectrum and p = (Q,u;)] to solve exactly the many- 
body equations. This elegant theory provides an explicit expression for the 
Green function valid only in the small energy range where the singularity 
dominates. This fact leads to problems in data analysis when the bare ex- 
citation energy u> (Q) reaches values well above 2A, since the signal around 
2A strongly decreases. To understand the correlation between the high en- 
ergy part of the spectrum and the one-excitation contribution, it is thus 
necessary to extend the validity of the theory to a wider range of energies in 
order to describe properly the continuum contribution to S(Q,uj). It then 
becomes crucial to consider the two main channels that contribute to the 
continuum: The excitation of a quasiparticle that then decays into a pair 
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of quasiparticles and the direct excitation of two quasiparticles. It is possi- 
ble to construct a phenomenological theory that takes into account this two 
channels. Details are given elsewhere,B here we only report the following 
simple expression for the density-density response function 



a 2 + 2a(3VzF{p) + fi 2 F{p) Q- l {p) 
X[P) ~ Go 1 ( P )-ViF( P ) 



In this expression (valid at zero temperature) a and (3 are the matrix element 
for the excitation of one and two quasiparticles, respectively. The coupling 
V3 parametrize the decay amplitude and F(p) is given by the sum of all the 
diagrams with two lines joined at the two external legs linked at least by two 
lines B The effect of a V4 interaction is completely hidden in F(p), and if we 
consider only energies near the threshold, Pitaevskii-ZRS theory can be used 
to obtain an explicit expression for F(p). The dynamical structure factor is 
simply related to x by S(Q,u) = — Imx(Q,w) and Eq. (|l]) constitutes the 
starting point to analyze data. We also note that Eq. ((!]) for /3 = gives the 
Green function of one quasiparticle. The above expression thus defines a Q 
and a x that share the same poles as expected in a Bose condensed systemS 



2. DATA ANALYSIS AND RESULTS 

Since the explicit calculation of F is a difficult task and in general 
depends strongly on the detailed structure of the vertex functions, we do 
not try to calculate it microscopically, but we extract it directly from data 
by exploiting the large (energy and momentum) region of validity of Eq. @. 
The main idea is to extract a one-dimensional function [-F(w)] from a two- 
dimensional experimental function [,S(Q,u;)]. As a matter of fact the main 
Q-dependence of Eq. (|l|) is through u> (Q), since F(Q,u) is expected to 
depend weakly on Q (as is the case for F (Q, to) in this momentum region). It 
is thus possible to extract F(u) and u; (Q) by fitting Eq. (|l|) to the different 
sets of data with different momentum at the same time. This procedure 
exploits fully the information contained in the data because it is sensitive to 
the correlation among sets with different Q. 

To extract the complex function F(uj) from data we parametrize its 
imaginary part with N real numbers (15 in the fit presented) in the following 
way. We choose a set of values of u, say {u>i, u>2, ■ ■ ■ , ujn} with 2 A = uj\ < 
u>2 < ■ ■ ■ < ujn reasonably spaced and we assign a free parameter, a^, to 
each of them. ImF(iv) can then be defined as a cubic spline interpolation 
on such a set. The real part can then be easily obtained by the relation 
ReF(uj) = —1/ttV J doj'Im.F (lo) / (u) — u)'). It is also important to reduce 
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the parameters to the minimum number of independent ones, so we define 
a dimensionless function f(uj) = XF(uj), constrained by the normalization 
condition fig tL)Tmf(u) = uj n - 2A, and g 3 = V3/A 1 / 2 , = j3/(\ l l 2 a). In 
this way a, j3, #3, and the N — 1 parameters that define / are independent 
and can be fitted to the data. 

We thus fitted Eq. (ffl) (convolved with the known instrumental reso- 
lution) to experimental data of Ref. 8 at 1.3 K. The resulting fit is shown 
in Fig. |]. The good agreement between theory and experiment is obtained 
with a reduced x 2 °f nearly 4, thus indicating that even if we are leaving a 
large freedom in / the agreement is significant. 




co (mcV) 



Fig. 1. Fit to the data from Ref. 8 with a parametrized Im/. In the inset 
the resulting Im/ is shown compared with ImF averaged over 2.3 < Q < 3.2 
A -1 and properly scaled. 

The new dispersion relation for the undamped excitation is reported in 
Table 1. Note that the value of 2 A is 1.484 meV.0 The fitted values of the 
normalized a% is also reported in Table 1. 

We find that the model can quantitatively explain that the peak position 
of S(Q,oj) is slightly larger than 2A for Q > 2.6 A -1 . This originates from 
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Table 1. Dispersion relation for the quasiparticle spectrum, weight (Z(Q)) 
of the pole, and energy of the bare pole u (left Table). Fitted values for a% 
at each uj; l (right Table). The values for Q = 3.6 A _1 are obtained from a fit 
not shown. 

a mixing within the instrumental resolution of the contribution of the sharp 
peak at energy ujq slightly smaller than 2A with that of the continuum of 
two rotons excitations starting at 2A. On general grounds the continuum 
should depend strongly on oj near 2A, as for co > 2A there are much more 
states available to decay into. Our procedure exploits the theoretical model 
for x(p) to extract both the value of ujq and the continuum contribution. In 
this way we can find the final part of the dispersion relation with improved 
accuracy. It turns out that data agree with a dispersion relation for the 
excitations always below 2A. 

Concerning the other parameters of the fit we find a 2 = 1.4, = —0.06 
meV -1 / 2 and (73 = 0.8 meV 1 / 2 . This implies that the direct excitations of 
two quasiparticles by the neutron gives a small but appreciable contribution 
to S(Q,lj), particularly at high energy. 

The shape of Imf found with the fit (see Fig. 1) has two main features: A 
clear peak at to ~ 2 meV and a "quasi-divergent" behavior at the threshold. 
The peak is due to the maxon-roton van Hove singularity as is clear by 
comparison with F Q . It is remarkable that although no trace of the peak is 
apparent in any of the experimental plots the procedure succeeded in finding 
correctly this information. 
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The quasi-singularity at threshold can be understood as an interaction 
effect, namely a signature of the roton-roton attractive interaction. As a 
matter of fact, in the small region of energy near the threshold we can apply 
Pitaevski-ZRS theory to evaluate F{oj). We verified quantitatively this fact 
by repeating the fit using Pitaevskii-ZRS theory to parametrize f(uj) and 
setting a cutoff in energy at 2A+0.2 meV. The resulting reduced x 2 f° r the 
fit gives evidence that the theory works quantitatively in this small region. 
We obtained in this way for the interaction parameter V4 « —4.7 meV A 3 
with a tiny bound-state energy of 1.3 jieV . Note that without the cutoff this 
theory is not able to fit data satisfactorily. 

In conclusion, we presented a theory for 5(Q,u;) that takes into account 
both one- and two-quasiparticle excitations by the neutron. The theory 
reproduces the experimental results over a large range of energy and mo- 
mentum. We have thus been able to extract the final part of the spectrum 
dispersion relation in 4 He and to determine its end as an hybridization with 
the bound state of two rotons. 
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